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I. INTRODUCTION 

A beautiful theorem by Leon Lichtcnstein demon- 
strates the existence of solutions describing the circular 
motion of two fluid bodies around their center of mass 
under the Newtonian force of gravity [l[. This solution 
is invariant under the symmetry generated by the vector 
field dt + fld<f,, the helical symmetry. 

Are there corresponding solutions in Einstein's the- 
ory of gravity? For truly isolated systems such solutions 
will not exist in General Relativity because the orbiting 
bodies will generate gravitational radiation and spiral to- 
wards each other. It is however conceivable — and this 
was first conjectured by Detweiler Q — that if there is 
the right amount of incoming radiation which would just 
compensate the outgoing radiation, then the bodies could 
stay on a circular orbit and the spacetime would admit 
helical symmetry as an isometry. 

There is an example of a solution with helical symme- 
try in Maxwell's theory by Schild [3|: He shows that two 
opposite electric point charges can move on a circular or- 
bit if the force on one particle is given by the Lorentz 
force of the "one half advanced plus retarded field" of 
the other particle. In this solution the total energy of 
the field is infinite. (See also [4| where this is generalized 
to N scalar point particles.) 

One also expects that in a fully relativistic treatment 
of the 2-body problem with helical symmetry the ADM 
mass of the solution will be infinite and that the solu- 
tions will not be asymptotically flat in the usual sense 
at spatial or null infinity. Because of this it is not even 
clear how one defines "helical symmetry" if there is no 
asymptotic symmetry group to characterize the infinites- 
imal helical Killing field dt + fic^ . 

Apart from this fascinating question whether Ein- 
stein's theory has really solutions with helical symme- 
try in which incoming radiation balances outgoing radia- 
tion, there is interest in the numerical community to use 



data with approximate helical symmetry in the numeri- 
cal construction of solutions of the 2-body problem and 
the calculation of the radiation emitted. This question 
has lead to various both analytical and numerical studies 
of helical solutions of linear and nonlinear model prob- 
lems by J. Friedman, R. Price, E. Gourgoulhon, their 
coworkers, and others (see, e.g., [|| - @, and references 
therein). The primary goal here is computing quasiequi- 
librium configurations of close binary compact objects, 
or individual rotating neutron stars. 

To address the problem of helical solutions in Ein- 
stein's theory it is absolutely necessary to have a really 
good understanding of the linear situation. The reason 
is that solutions to nonlinear problems are usually solved 
as a limit of a sequence of linear problems. No system- 
atic study of linear helical solutions of the scalar wave 
equation, Maxwell's equations and linearized gravity is 
available. This is the topic of the present paper. 

Section II studies solutions with helical symmetry - 
i.e. solutions invariant under the symmetry generated 
by dt + Qd<j, - of the scalar wave equation on Minkowski 
space. There is a very important point: because of the 
assumed symmetry one can study the wave equation on 
the space of orbits of the symmetry. There one obtains an 
apparently simpler 3-dimensional problem. However, the 
symmetry reduced wave equation is a second order PDE 
which changes its character: it is elliptic near the center 
and hyperbolic far out. No general theory is available in 
this situation. There is the paper by Torre [lfj| who ob- 
tains some existence results for the inhomogeneous equa- 
tions on part of Minkowski space. 

We proceed completely differently. A lot is known 
about the wave equation in Minkowski space. In partic- 
ular, the retarded and advanced solutions of a spatially 
compact source with helical symmetry are also helical. 
So, with almost no work we obtain solutions of the inho- 
mogeneous wave equation on the space of orbits. 
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The general homogeneous solution with helical sym- 
metry can be constructed via a spherical harmonics de- 
composition. So, we know the general solution. 

There are helical solutions with explicit Yi m angular 
behaviour. We give an integral representation of such 
solutions. 

In Section III we investigate the asymptotic behaviour 
of these solutions near spatial and null infinity. The re- 
tarded solution of a spatially com pac t helical source has 
the usual Bondi-type expansion at future null in- 
finity but no such expansion exists at past null infinity. 
The standard radiation field of the solution, as defined, 
e.g., in [ll[, pjj], [H[> does not exist on past null infinity. 
Hence for the retarded plus advance solution no radiation 
field is defined at all. Also, the decay at spatial infinity 
is rather weak, the field decays like sin r/r. 
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(In Sees. II and III we set the speed of light c = 
1.) A remarkable property of this equation is that it 
changes its character! The norm of the helical Killing 
vector is N = — 1 + Q 2 r 2 sin 2 9 and therefore the he- 
lical Killing vector is timelike near the axis (N < 0), 
null for N = — 1 + f2 2 r 2 sin 2 6 = 0, and spacelike for 
N = -1 + tt 2 r 2 sin 2 8 > 0. This implies that the re- 
duced equation ((2|) is elliptic near the center and hyper- 
bolic for N = — 1 + 2 r 2 sin 2 9 > 0. This raises interest- 
ing questions concerning solutions defined globally. We 
will denote the three regions as: T = {(</), r,0)\N < 0}, 
L = {($, r, 8)\N = 0},S = {(0, r, 8)\N > 0}. 



In Section IV we calculate the retarded and advanced 
field of a scalar source moving on a circle with constant 
velocity. The rather complicated periodicity properties 
of the solutions become apparent. 

The analogous solution in Maxwell's theory is pre- 
sented in Section V. We give the fields in terms of the 
null tetrad NP (Newman-Penrose) projections, i.e., as 
three complex scalars. Since the charged particles mov- 
ing in circular obits occur in synchrotrons (where they 
are guided by external magnetic fields) our final expres- 
sions might be of some use also elswhere. In Schild's 
work 0] the field is not given. The usual radiation field 
again does not exist but the NP scalars exhibit general- 
ized peeling properties. As a consequence of the helical 
symmetry, the retarded and advanced NP scalars show 
specific symmetry properties. 

From the results in Maxwell's theory one can almost 
guess the formulae for Einstein's linearized gravity if 
needed. 



II. THE WAVE EQUATION ON MINKOWSKI 
SPACE 



We begin by considering the wave equation on 
Minkowski space and want to describe global and local 
solutions with helical symmetry. 

A function F(t, r, 9, </>) (in standard spherical coordi- 
nates on Minkowski space) is invariant under dt + Qd<p if 
and only if 



F(t, r, 9, </>) = G{tf> - nt, r, 9) = G($, r, 9) 



(1) 



The function G is periodic in its first argument with pe- 
riod 2it. F is a solution of the scalar wave equation if and 



A. Local source— free solutions 

Obviously in T local solutions can be determined by 
boundary value problems and in S by Cauchy problems. 
What happens near L is not so clear. 

Suppose we have a C 2 solution of ^ defined for n < 
r < T2 such that N changes its sign in this region. We 
can expand into spherical harmonics 

G(lr,6) = G(r)Y lm 0,6) , 0e[O,27r), (3) 
and obtain for G(r) the radial equation 



m 2 n 2 G+ ~{rG) 
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T l(l + 1)G = 
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The general solution in terms of spherical Bessel func- 
tions is 



G lm {r) 



" hn -.J l+ >(rnnr) + J^=Y l+ >( m nr). (5) 



The set L, the light cylinder, is no singularity for the ra- 
dial equation. Hence, we see that there are local solutions 
for which nothing particular happens at L. 

The only solutions regular at the origin are given by 
b = 0. For mf2r 3> I the solutions behave as sin(mf2r — 
fZ)/mOr. 

Thus we know the general C 2 solutions in terms of con- 
verging superpositions of spherical harmonics and Bessel 
functions. 



G(4>,r,9) 
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The coefficients a; m , 6/ m in {5| may depend on f2. Taking 
([5]) in the form 



G lm (r) 
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where a; m ,6/ m are independent of f2, we find, using the 
properties of the Bessel functions, that in the limit — > 



Gi m {r) 
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(8) 



This is the general static Yj m -solution of the Laplace 
equation. 



B. Global source— free solutions 

Now we assume that we have a C 2 solution defined on 
Minkowski space. Again we can decompose into spherical 
harmonics and because of the regularity at the origin only 
J-Bessel functions appear: 
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This shows that a global C 2 solution is uniquely deter- 
mined by its values in a neighborhood of the center. 

What is the asymptotics of such solutions? From the 
Yim decomposition it is not obvious what type of decay 
one can produce by superpositions. Remember that a 
superposition of not decaying plane wave solutions of the 
wave equation can produce solutions with spatially com- 
pact support. 

There is an easy argument that we cannot obtain so- 
lutions of finite energy: suppose we have a C 2 solution 
of finite energy; then the energy in a finite box around 
some fixed point in space has to decay in time, a general 
property of solutions of the wave equation. Such a decay 
contradicts helical symmetry. 

In particular, there are no C 2 solutions with spatially 
compact support and helical symmetry. 

In the next section we shall see that there are local 
solutions near null infinity which allow a Bondi expan- 
sion, but global source-free solutions have never a Bondi 
expansion. 



C. Global solutions with sources 

Suppose we have a solution of the form ([9|) defined near 
the center. In general, the series will not converge for all 
r. Examples are given by solutions of the inhomogeneous 
wave equation with sources invariant under the helical 
symmetry. The simplest case is a point source moving 
on an integral curve of the helical Killing field. In space- 
time the retarded and advanced fields of this source are 
well defined outside the source. Near the center these so- 
lutions must be of the form ([9]) but they become singular 
at the source position. Further out we will again have a 
superposition of now both Bessel function solutions. 

The general solution of the inhomogeneous wave equa- 
tion with a spatially compact source is given by the sum 
of the retarded solution and some global solution of the 
homogeneous equation. 



If the spatially compact source has helical symmetry 
the (unique) retarded and advanced solutions have also 
helical symmetry. 

A simple way to see this is to consider a point charge 
as a source, as we do in detail in Section IV, and change 
the "source variables" t' , (/>' to t' + A and 4>' + £1A, A G R, 
which represents the action of the helical symmetry on 
the source. It is then easy to show, using the explicit 
formulae for Lienard-Wiechert potentials, that the field 
remains the same when the "observation variables" t, <p 
change to t + A, tj> + £1A. For an extended source the 
superposition principle can be used. There is, in fact, 
a more general theorem available showing for the wave 
equation that a symmetry of the source is inherited by 
the retarded/advanced fields. 

If we take these helical solutions, which we obtain from 
general theorems about the wave equation in Minkowski 
space, and consider them as solutions of the reduced 
equation (2), we obtain immediately an existence the- 
orem for the reduced equation which, in contrast to the 
result by Torre (l0| . includes the center. 

More precisely: The retarded and advanced solution of 
the wave equation with a spatially compact helical source 
defines a global solution of the inhomogeneous reduced 
equation which is analytic outside the source. The gen- 
eral solution is obtained by adding a global solution of 
the homogeneous equation described in Section II B. 

If we take a source with helical symmetry, then the 
difference of the retarded and advanced solution defines 
a global, source-free solution with helical symmetry on 
Minkowski space. This shows the existence of source-free 
solutions defined on all of Minkowski space. 

We want to give an explicit representation of a helical 
solution with Yi m -behaviour: Let p be a helically sym- 
metric source of the scalar field r), satisfying the 
inhomogeneous wave equation 

□ $ = -4wp{t,r) , (10) 

where r = (x l ) = (x, y, z). The retarded solution reads 



r-r 



with 



t = t- |i 



— r| = y r 2 — 2rr 



(11) 



(12) 



(13) 



r = |r| = (x l x l )2 , similarly for r. 

We assume for the helical source the Yi m dependence: 

p(4>,r,6) = p{r)Y lm {4>,6) = (14) 

= p{r)ai m e" n *Pr (cos 9) , $ G [0, 2ir) . 

The integrand / of the retarded integral becomes 

/ = r -^ =r p(f)ai m e im [*- n (*-l r -*IMiY B (oose) . (15) 
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So, 



A. The scalar wave equation 



p(r) 



I = -^-^-a im e jm<#, e-™ n *e jmn|r - F| P i m (cos^) , (16) 



/ = e - mnt >(f)y Im (0, 



tOjr— r| 



r - r 



(17) 



Now we can use the following expansion in spherical har- 
monics (see [lj], Eqs.(16.18)-(16.22) with mO = k): 



4tt\t — r| 



= J2m(r,f)Y^(S,m m (0,(l>) , (18) 



I, in 



g(r,r) = ikji{kr < )hi { - r> {kr > ) , 



(19) 



with r < — min(r,r),r > — max{r,r). The Bessel and 
Hankel functions ji, hi^ can be found in [l4| [Eqs.(16.9), 
(16.10)] or [H]. 

If we insert this expansion into the retarded integral, 
use the orthogonality properties of the spherical harmon- 
ics and remember the volume element we obtain 



* = e- mQt Y lm ( 



(20) 



x / drr 2 p(r)imClji(mQ,r < ^)hi^ (mf2r>) , 



$ = e-"" ot r i „ l ((/),6')imfi/i; (1) (mr2r) dff 2 p(f)ji(mfif) 



e- imnt Yi m {(j),0)imnji{mnr) \ dff 2 p{f)hi {1] '(mOf 



(21) 

For compactly supported p(r) the behaviour near infinity 
is given by /^^(raOr) ~ e %mQ " r jmSlr. In this represen- 
tation of the solution the "retardation in space time" is 
completely hidden. 

Thus, we have complete control of all solutions of the 
inhomogeneous reduced equation with signature change 
on the quotient. The lesson is not to forget Minkowski 
space and just work with the reduced equation on the 
quotient of the symmetry, but use all the results one has 
on the wave equation in Minkowski space! 



III. ASYMPTOTICS OF SOLUTIONS WITH 
SOURCES 

In this section we discuss the asymptotics of helical so- 
lutions on Minkowski space for scalar fields, and we also 
mention the analogous cases for Maxwell fields and lin- 
earized gravity. Here we deal with general (unspecified) 
sources of spatially compact support. In the next sec- 
tions we investigate the fields of point sources on circular 
orbits in detail. 



Let p be a helically symmetric source of the scalar field 
$(t, r), satisfying the inhomogeneous wave equation 



□ $ = -Anp(t,r) 



(22) 



where r = (x l ) — (x,y,z). Again, the retarded solution 
reads (fTTj) . Assume now general p which is bounded in 
space; then at large distance from the matter, r > f, we 
obtain expansions for p in (fTTj) . 
We write 

\T-T\ = r + a(f ) <f>,eJ ) S) + o(^\ , (23) 



and 



r — r r \r 



Because p is helical we have 

p(i,$,6,f) = p$-Sli,6,f) , 



(24) 



(25) 



where p is periodic with period 2ir in its first argument. 
Hence we obtain 

4M , M= /fc3^t3iM^. (26) 

J l r - r l 

Let us first investigate the field for r — > oo, t fixed. We 
can expand, using (f2"3l) . 

p$-n(t-\v-T\),9,f)= (27) 

= p$-n(t-r-a(f,<p,0,$,6))J,f) 

l|(^-fl(i-r + a(f,^ ) ^l) ) S»xoW+... 

(We can easily write a complete Taylor expansion if 
needed.) 

The leading contribution to the field is therefore [see 
also ([Mjl ] 



$(t,<l>,6,r 
1 



(28) 



p (4> - Cl(t - r - a(f, cj), 6, </>, 9)), 6, f)dV+0 



The function 



P(t-r, <f>,0) = I p $ ~ 0(* - r - a(f, </>, 6, & 9)), 0, f) dV 

is periodic in t — r with period 2tt/CI and the asymptotic 
behaviour is oscillatory in r for t fixed, r — > oo: 



9>{t,<l>,6,r) 



P(t - r, , 



(30) 
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For a helical source with Yi m behaviour, 
p($,6,r) = e im *P l (cosd)f(r) 



(31) 



we can see the periodicity explicitly. 

Next we consider the asymptotic behaviour at future 
null infinity T + . With u = t — r, we have 

p($-tt(u + r-\r-r\),6,r) = (32) 
= p(0-fi(u-a(f,&M,0)),0,f) + 

+|(^-n(ii-o(F ) ^M,j)) 1 «»x0^j+... . 

If we put this in the retarded integral ([TTjl we find that 
for u =const, r — ► oo, we obtain the usual Bondi-type 
expansion in r _1 (see [ll|, [TH) as we approach X + . 

If we go to past null infinity T~ , we have v = t + r and 
obtain 

p (4> - Q(v + 2r - \r-f\),e,f) = (33) 
= p (4> - Q(v + 2r - o(r, <£, 0, & §)),§, f) + 

+|| - n(« + 2r - o(f, 0, 6, &)),§, r)xo(~^+.... 

This for r — > oo, u =const implies again oscillations. This 
behaviour is easy to understand in a spacetime picture: 
when going to Z~ through the outgoing field of a source 
which has been periodically moving at all times, we cross 
infinitely many oscillations. 



B. The Maxwell field 

We take a continuous compact distribution of charges 
with helical motion. They define a conserved, helical 4- 
current. If we use Lorentz gauge we have for each com- 
ponent of the 4-potential to solve a scalar wave equation 
with a helical source. Because the source is spatially com- 
pact the Lorentz condition is satisfied for the solutions 
given by the retarded, or advanced potentials. Hence, all 
the results for the scalar field apply. 



C. Linearized gravity 

If we take the T^ v of particles moving along the orbits 
of a helical Killing vector we can consider the linearized 
equations in harmonic gauge. We then encounter again 
wave equations for individual components and can write 
down retarded and advanced solutions. The asymptotics 
is thus very similar to that of the scalar fields. However, 
now we cannot satisfy the harmonicity condition because 
the energy momentum tensor in linearized gravity for cir- 
cular orbits is not conserved. A difference to the Maxwell 
case. 

If we want a complete solution the matter must satisfy 
the equation of motion. For example, those of linearized 
elasticity. 



IV. THE FIELD OF A POINT SCALAR 
CHARGE IN A CIRCULAR ORBIT 

We adopt the standard procedure of finding retarded 
and advanced Lienard-Wiechert potentials in electrody- 
namics ( see, e.g., an elegant, covariant description by 
Rohrlich [TfJ]) to the case of scalar fields. The retarded 
(e = +1) and advanced (e = —1) solutions of the wave 
equation with a 5— function type source with scalar 
charge Q moving along the worldline x**(t) so that 



p(x") = Q I S (4 \x» - x^(t))oIt , 



can be written in the form 



.9. 



where 



w, 



-e Va/3 



(34) 



(35) 



(36) 



(Henceforth, we do not set set c = 1.) Here the null 
vectors 



(37) 



connect a given spacetime point x@ , in which $ c is to be 
calculated, with two points if on the worldline x@(t), in 
which the past and future light cones with vertex at x 13 
intersect the worldline. The 4-velocities vf are evaluated 
at these points if. The null vectors can be written as 

= (eE = e|R|,R=r-r e ) , (38) 

the time components R® — eR are equal to the retarded 
and minus advanced distance of the particle to the field 
point x a . 

Now consider a point charge Q moving in the plane 
z = along the circular orbit of radius a with center at 
x — y — z = 0. So 



a (cos(fit + 0o ), sin (fit + O ), 0) 



(39) 



where f2=constant is the particle's angular velocity, 
is the azimuth at which the charge occurs at t — 0. The 
4-velocity is 



= (c7, tv) 
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1 



a 2 Q- 



constant , 



(40) 



(41) 



the 3-velocity is 

v = afi(-sin(Qt + </>o),cos(fit-|-0o),O) . (42) 

From these simple relations we find the distance R en- 
tering the null vectors (|3"T|) to be given by 



Re 



( (x — acos(ftt e + </>o)) 2 

+ (y - a sm(nt e + <?(>o)) 2 + z 2 ) 5 



(43) 
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which, in spherical coordinates, implies 



R e = r ( 1 — — sin#cos(</) — f2f e — <po) + —. 

r r 



2 \ 2 



(44) 



Retarded and advanced times are implicitly given by 
the equation 



t f = t-e- 



(45) 



V ( 2ft (X 

t — e- 1 sin 9 cos(0 — flt e — d>o) H — 

c \ r r 



2\ 2 



Note that on geometrical grounds we know that this 
equation has a unique solution for t t provided the world- 
line of the source is timelike, i.e., — < 1. In contrast 
to the case of linear uniform or of uniformly accelerated 
motion, it is however impossible to obtain an explicit ex- 
pression for t e as a function of (t, r, 9, </>). 

The expression in (|36|) for the scalar u> £ now reads 



w e — jr [ ( 1 — — sin#cos(0 — flt e — O ) + ° 



2\ 2 



-e — sin0sin(0-ftt e - cf> Q ) ] , (46) 



(53) 



For a subluminal source, a < 1; in the non-relativistic 
limit, a<l. At large r»a we find 



^ = 1-^ + 1^^+011- 



(54) 



irrespective of t,t e ,6,cp. This implies for the scalar field 
to be given at r 3> a by the expansion 



$ f = 



77- (1 — ea<S e 



1 



Ce a 



(1 — ea5 e 



(55) 

For = we of course get the asymptotic static field of 
a source at rest at r = a, 4> — 4>o, = 7r/2. 

For a moving source we have to take into account equa- 
tion (|45[) for the retarded/advanced time by employing 
the asymptotic form ([54)) for TZ e : 



t E = i - e-U c =t-e- + e-C e - e— S 2 — + O — 
c c c zc r z \r z 

(If we drop the O-term we obtain a new, approximate 
implicit equation for t e .) Fixing t, 9, 4>, we get 



which then implies directly the form of the resulting 
scalar field ([55)1 . 

Let us finally rewrite the last relations in terms of the 
'corotating ' angular coordinate <f> = <f>~ fit. Introducing 



(47) 



the equation for the retarded/advanced time becomes 



fir / 2a a 2 \ 2 

> + e — 1 sin0cos((^- <M + -( 48 ) 

c \ r r z ' 



Then the resulting scalar field, 



Ol/, 2a . „ , N . a 2 ^ 2 



ill sin6>cos(0 e - O ) + ., 

7r I y ?" 7* 



afi 

-e — sin0sin(0 e — </> ) 
c 



(49) 



is independent of the time t. 

Let us finally investigate the asymptotic behaviour of 
this solution. Define the following dimcnsionlcss quanti- 
ties 



C e = sin0cos(0 — f2i e — 00 ) j 



iS e = sin#sin(0 — fi£ f 



2 \ 1/2 

' R - = I 1 " 2C ^ + 72 



(50) 



(51) 



(52) 



<9i e 
<9r 



-e/c 



t,8,<j> 



1 



ea5 f ' \ r 



(57) 



Since 1 — eaiS e > 0, t e is monotonically decreasing (e = 1) 
or increasing (e = —1) function of r. 

Consider first the limit t fixed, r — > 00: The functions 
C e , S e are oscillatory as r — > 00 and, correspondingly, the 
scalar field decays, while oscillating, in accordance with 
Eq. ([53)) . Retarded and advanced fields have similar 
asymptotic behaviour at spatial infinity io. 

To see the character of the oscillations in more detail, 
we solve the equation for the dimensionless quantity flt e 
iteratively (realizing that the term |eaC e | < 1). So, 

flt t ~ flt-e— + easing cos(0- 6 - fi(i-e-)) . (58) 
c c 

The functions appearing in the fall-off of $ c , with fixed 
t, 9, <f> and r — > 00, thus behave as 

f n . n \ 

St — pi sin I e— r + ep 2 cos(e— r + P3) + Pi J , (59) 

Pi . . - Pi are constant parameters depending on fixed 
t,6,(f>; similarly for C e . Notice that S e is again periodic 
in r with the period 2nc/fl. 

At null infinity, retarded and advanced fields behave 
very differently. Turning to future null infinity T + , we 
put u = t — - and consider the limit r — > 00, with it, 9, (j) 
fixed. From the asymptotic form of the equation (|56() for 
t f for r — > 00, u, 9, <fi fixed we now find that -+- = 0(r _1 ), 



whereas 



dt_ 

dr 



: (1 — aS- ) 1 . Here ± is the abbreviation 
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for e = ±, i.e., retarded/advanced. Hence, t+ is approx- 
imately constant at large r, but i_ is a monotonically 
increasing function of r. Indeed, using the retarded time 
t + from ([55]) , we get C+,6>+ as functions of u,6,(f) , so 
that the retarded field given by (|55[) with e = 1 has at 
X + the standard Bondi-type expansion [ll|, [l2j 



$, = — 



7 



h{u,e,4>) , / 2 (M» 



1 



However, the advanced field — 1[55]) with e 
contains functions like 

/ 2 2 
<S_ = fei sin I — -f2r + &2 cos( — fir + fca) + k± 



(60) 
-1— 

(61) 



which oscillate in r with period 7rc/f2, so half of that at 
spatial infinity. 

At past null infinity I~, r — > oo with i? = t + ^-,0^ 
fixed, the advanced fields exhibit the standard Bondi 
type expansion whereas retarded fields decay oscillating 
in r. (All this becomes very clear in a spacetime picture.) 

The oscillating factors at r~ k in the asymptotic form 
([55| of <3?e can become large for ultrarelativistic velocities 
when the source moves near the light cylinder. 

For non-relativistic velocities, the asymptotics of the 
fields simplifies considerably. Assuming a <C 1 and ne- 
glecting terms which are 0(a 2 ), we find the asymptotic 
expansion ((55)) to yield 



— |l + easin^sin (j) — 0o — — e_ )) } 



--^ |l + e2asin6»sin (j> - (f> - fl(t - e-)) } 
x {sin 6> cos 

From this simple explicit expression the asymptotic be- 
haviour of 4> e at io, 1^ becomes transparent. In partic- 
ular it is easy to see that at the leading order near ?o 
both $ e oscillate as r _1 sin(efir/c), whereas at I + the 
advanced field behaves like r™ 1 sin(2f2r/c) and, in the 
same way, the retarded fields fall off at T~ . 



V. THE FIELD OF A ELECTRIC POINT 
CHARGE IN A CIRCULAR ORBIT 

As mentioned above, in the Lorentz gauge the results 
for the individual components of the 4-potential will be 
essentially the same as for the scalar field. It is, however, 
of interest to see the behaviour of the electromagnetic 
field tensor because it determines directly the force and, 
from the perspective of the asymptotics, the peeling- 
off properties of its null tetrad projections. Since these 
may find applications also in other context, we shall give 
the null-tetrad components in the whole spacetime for a 
point charge moving along a circular orbit with, in gen- 
eral, relativistic speed (and so emitting synchrotron ra- 
diation). 



We start from a general, covariant expression for the 
retarded/ advanced field tensor of an electric point charge 
e moving with a 4-velocity v a and a 4-acceleration a a 
(l6| . It is convenient to define a unit spacelike vector u a 
which is orthogonal to v a , u a u a — 1, u a v a = 0, such 
that the null vectors R% = x a - x% (see Eqs. J37]), (|5g])l 
can be written as 



Conversely, 



(63) 



(64) 



where the last expression is identical with (l36|) . The 4- 
potential is simply 



C W e 



(65) 



which shows that it differs from the scalar field <& e given 
in ([33)1 just by the replacement Q — ► — ~v a . The resulting 
expression for the Maxwell field tensor reads: 



e 



(66) 



where the small square brackets denote antisymmetriza- 



tion without " |" , a e 



For the circular motion with the same parameters as 
with the scalar charge all expressions for the velocity, 
retarded and advanced times and the scalar w e are again 
given by formulas (|3T))) - (|4"B"|) . In the expression (f6"6"|) we 
also need the 4-acceleration a M = dv^/dr = 'ydv^/dt, 
where v 11 is given by (|4"0"]) - (|42)k the result is 

an 2 

(O,cos(m e + o ),sin(m e + (/>o),O) . 



1 



a 2 Q 2 



(67) 

From (|6"3"|) the spacelike vector u£ can be expressed as 
u% = R^w e ~ x — ev^/c, which in case of the circular mo- 
tion implies the components to be 



u: — eu> 1 rR f 



(68) 



w e 1 (r sin cos — acos(flt e + 0o)) + (69) 

+ -7<if2sin(f2t e + 4>q)) , 
c 

w^ 1 (r smOs'mtp — asm(flt e + 0o)) — (70) 



-~faQ cos(Sli c + tfio)) 



w c r cos f 



(71) 



it is easy to check that is unit and perpendicular to 
v£. The last necessary ingredient is the projection a e of 



ar on u 1 ^ 



-wrVftV sir 



(72) 
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The 4-potential is determined by (|65j) in terms of the 
quantities we know already from the scalar field case — 
the scalar w e , see (|46]l . and by the 4-velocity v£ given by 
(f4"0")) - l[4"2")) with t — ► t e . Notice that the time-component 
A\ is equal to the scalar field <J? e [see ([35])] up to the sign 
and the multiplicative 7 factor. 

For highly relativistic motions, A\ is thus much big- 
ger than the corresponding <& e . The asymptotics of A\ 
is, however, the same as of $ e (see (|55|) . (f60|) and (|62|) V 
Owing to the symmetry, A\ — 0. Because of the velocity 
entering the componets A^ , A^, their asymptotic expan- 
sion at io begins immediately with the oscillatory terms 
at r _1 even for non-relativistic motions (a <C 1), rather 
than with constant xr _1 as in (|62j) . and the same with 
the advanced (retarded) fields at I + (I~ ) . 

Substituting for w e ,v^,a^,u^ into dHS]) we obtain all 
components of the Maxwell field tensor. These are quite 
lengthy, here we give just the null tetrad components 
from which all information can be retrieved. The stan- 
dard null (Newman-Penrose) tetrad (see, e.g., [17[) in 
Minkowski spacetime reads |18| 

I* 1 = — cos 0sin0, sin sin#, cos 9) , 
v2 

= — — cos0sin#, — sin0sin#, — cos#) , 
v2 

m p = -^(O,cos0cos# — isin0, (73) 

sin cos 9 + i cos 0, — sin 6*) ; 

in our signature it satisfies l^n^ = — 1, m M m M = 1, all 
vectors are null. 

Now the null tetrad electromagnetic scalars are given 



as the following projections of the Maxwell tensor: 



$0 
$1 

$2 



1 

2 



(74) 



The null tetrad naturally induces the orthonormal tetrad 
(see Ref. [18]) associated with the spherical coordinates 
with the standard orthonormal triad ef,eg,e^ (supple- 
mented by the timelike vector (1,0,0,0)). From the 
scalars $ 9 , q = 0, 1, 2 the physical components of the elec- 
tric and magnetic field in the triad can be expressed by 
means of the relation 



Ef — iBf. = 2(f>i , 

E § - iB § = -$ + $2 , 

E A - iB i = -z(<t> + $ a ) . 



(75) 



The Maxwell tensor field is given by (|66|) . with all the 
ingredients expressed explicitly in (j3"9")) - (|l6")) . ([ST]) - ([72")) . 
In addition to C e ,S e ,TZ e and a determined by ([50]) - ([53]) . 
we introduce another dimensionless quantity 



Q e = lZ e — eaS e = 



„ a a 
l-2C £ - + -2 



eaS f 



(76) 



Finally, substituting the explicit expression for F^ u and 
the null tetrad (|75|) into (j?4"|) , we obtain -after straight- 
forward though not very short calculations and arrange- 
ments - the following results for the 3 complex scalars 
<&o, $1, $2, in which again C — — f2£ £ : 



$1 

$ 2 



2crQ 3 f 



„ „ o. or 
-a I 1 - eK e - C e - + -= 



cost/ cos 



0o ) - isin(0 £ - 0o ) 



(77) 



00 ) - JCOs(0 e - 0o ) 



a^(l - eft e -C e -) 



a. a 



+ 



2(rQ £ ) 3 
e 



COS ft COS( 0e — 00 ) — « Sm( e — 00 



2r 2Q3 

ef2 



1 — ea7?. f sin 9 sin( 



2crQ3 



-a(l + eTZ e -C £ - + - 



- ( crC e + (1 - a ) - + ia cos6>) sin#cos( 



cos V cos 



00 ) + «sin(0e - 0o ) 



+ (l + eTZ e )- 



2(rQ e )3 



COS V COS I 



- O )- ZCOS(0 £ - 0o ) 

bo) - isin(0 e - 0o ) 



a 2 (l + e^e-C e -) + - 



0o ) +ia-cos6>) ,(78) 
r ) 

(79) 



(80) 
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There are well-known properties of the null-tetrad com- 
ponents of the zero-rest mass fields (see, e.g., In 
the electromagnetic case $ 2 characterizes the outgoing 
radiation field at 2+: it is given as lim r ^ 00iU=co „ st (r<f>2). 
At T~ lim r _ >00iV=con , s t(r$o) determines the incoming ra- 
diation field. We shall see that, indeed, this is the case 
for $2+ at I + and for $o- at I~ . (Here ± is again the 
abbreviation for e = ±1, i.e., retarded/advanced fields.) 

Before turning to the asymptotic properties of the <&s 
let us observe some intriguing "symmetry relations" of 
these quantities. They arise as a consequence of the he- 
lical symmetry of the source and the fields. In order 
to understand them we have to turn back to the equa- 
tion (|45l) for retarded and advanced times. Notice first 
that putting the 'field' time t — 0, Eq. (j45]) still cannot 
be solved explicitly for t e , however it is easy to see the 
following relations between t + and t- : 



t+(t = 0, r, 8, 0) = -t- (t = 0, r, 6, -0 + 20 o ) 



11) 



(It is instructive to draw the circular orbit of the source 
and convince oneself that this relation can also be under- 
stood on "geometrical grounds".) Owing to the helical 
symmetry the relation (|81|) can be generalized for any t 
to read 

t+(t, r, 8, 4>+m)-t = -t-(t, r, 8, -0+f2i+20 o )+t ; (82) 

this can be checked directly using |45|) again. 

Now as a consequence of helical symmetry and, in par- 
ticular, of (|8"2")l . there arise the following relations — valid 
at any spacetime point — between the retarded and ad- 
vanced fields as given by the null-tetrad components of 
the $'s: 



$o+(t,r,6,4> + SH) = -$ 2 - 0,-0 + fit + 2<fo) 



<f> 1+ {t,r,6,<p + Q,t) = +$i_(t,r,0,-0 + Of + 20 o ) . 

It is straightforward to convince oneself that our result- 
ing expressions (j77|) - (f79|) satisfy these relations. They 
are the keystone in Schild's 2-body problem, in fact the 
iV-body as well, implying that tangential forces on parti- 
cles vanish when both retarded and advanced effects are 
taken into account. But this problem will be considered 
elsewhere (for the ./V-particle problem in the case of the 
scalar field, see Q). 

We now turn to the asymptotic properties of the fields. 
As with the scalar field we first make the expansions of 
the retarded and advanced fields at large r»a using the 
expansion (|54p of 1Z t . In addition to the the abbrevia- 
tions C e ,S e defined in (1501). (51]), i.e., 



C e = sin0cos(0 e — 0o) , (84 
S e = sin 8 sin(0 e — 0o) , (85 

we introduce 



C t = cos6*cos(0 e — 0o ) , (86 
S e = cos#sin(0 e — 0o ) , (87 

and 

C e = COs(0 e - 0o ) , 

s e = sin(0 e - 0o ) . 



$ 2+ (f,r,0,0+fii) = -$ O -(*,r,0,-0+m + 20o) , (83) 



With these notations we obtain the following: 



$2+ 
$0+ 



;(1 - aS+) 3 



+ is + — ia sin Q^j — hO 
{C_ +is- + a[C_(<S_ -icJ] 



2(1 + aS-) 3 

-i sin 8 - ~a(l +C 2 _)(C- +is) -C 2 _) sin 6>]}4r + (\ 

2 2 r 6 \r 4 

e 1 / 1 

[1 - eaS e - (a 2 C e + ia 2 cos 9)C e ] — + O — 



2(1 - eaS,) 



:{C + — is+ + a[— C+(S+ - ic-i 



2(1 - aS^ 

+isin8 - ~a{l + Cl)(C + - is+) - iy (1 - C 2 + ) sin0]}i + O 



$o- = - 



efla 



c(l + a<S_) 3 



C_ — is_ — ia sin 8) — \- O ( — - 



(90) 

(91) 
(92) 

(93) 
(94) 



10 



These formulas get still simplified for nonrelativistic 
motion of the source, a« 1, when only terms linear in 
a are kept (using (1 — eaS e )~^ ~ 1 + 3aS e + 0(a 2 )) and 
rearranging the terms: 



$ 2 + = ( C A 



(95) 



$2- 



ea 

Y 



(1 - 2aS-)(C- + is-] 



$ie = |(l + 2ea5 6 )i + o(i 



(96) 
(97) 



<J>o + = y 



(1 + 205+)^ 



3 + °fi) ' < ! »' 



$n_ = 



eOa /- \ 1 

C- — ?S— 



•O 



(99) 



functions of u, 9, <p [see Eq. (|60|) and the text therein]. In 
addition, in the electromagnetic case we observe a new 
feature in the asymptotics - the well-known " peeling-off 



property" [12J, [17J: $ 9+ ~ /,(u, 6>, </>)r?- 3 . The limits 



r 3_9 $ 9+ are well defined at r — + oo with u = £ — r/c, w, 
fixed, i.e., at With the advanced field at I + the sit- 
uation is different. We still find a "generalized peeling" 
in the sense that <E> 9 _ ~ r~ 1-9 ; so $ 9 _ corresponds to 
^(2-?)+- However, the limits r 1+9 $ 9 _ at r — > oo,u, 
fixed do not exist since these are oscillatory functions as 
<S_ in (|61[) . i.e., with period nc/Q. which is half of that 
at spatial infinity. Nevertheless, notice that r 1+9 <E> 9 _ re- 
mains bounded as r — ► oo. 



As expected, at past null infinity T~ , r — > oo, u = i + 
r/c,6,(j> fixed, the advanced fields exhibit the standard 
Bondi-type expansion and peeling, whereas the retarded 
fields do decay with r — > oo but in an oscillatory manner. 

All these features become evident especially when re- 
garding the simple fields (|95[) — (f99|) in the non-relativistic 
limit a«l. 



Clearly, the fields have again the oscillatory character. 
The retarded and advanced times t e can be expressed in 
terms of the field time t and functions C e and S € according 
to Eq. (f56|) . At fixed t,9,4> and r — > oo, i.e., at io, the 
oscillations are given by C e , S e which have the form |59|) 
as in the scalar-field case. The retarded and advanced 
fields behave similarly. 

At future null infinity I + the retarded field has the 
standard Bondi-type expansion since C+, 5+ become 
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